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Small black holes in heterotic string theory have vanishing horizon area at the supergravity level, 
but the horizon is stretched to the finite radius AdS2 x S D ~ 2 geometry once higher curvature 
corrections are turned on. This has been demonstrated to give good agreement with microscopic 
entropy counting. Previous considerations, however, were based on the classical local solutions 
valid only in the vicinity of the event horizon. Here we address the question of global existence of 
extremal black holes in the D-dimensional Einstein-Maxwell-Dilaton theory with the Gauss-Bonnet 
term introducing a variable dilaton coupling a as a parameter. We show that asymptotically flat 
black holes exist only in a bounded region of the dilaton couplings < a < a CI where a CI depends 
on D. For D > 5 (but not for D = 4) the allowed range of a includes the heterotic string values. 
For a > a CI numerical solutions meet weak naked singularities at finite radii r = r cusp (spherical 
cusps), where the scalar curvature diverges as \r — r cusp | -1 / 2 . For D > 7 cusps are met in pairs, so 
that solutions can be formally extended to asymptotically flat infinity choosing a suitable integration 
variable. We show, however, that radial geodesies cannot be continued through the cusp singularities, 
so such a continuation is unphysical. 

PACS numbers: 04.20.Jb, 04.65. +e, 98.80.-k 



I. INTRODUCTION 



During recent years important progress has been achieved in understanding the entropy of the so-called small black 
holes (for reviews see [USUI]) which have vanishing horizon area at the supergravity level [1, [1,(1,0] • The discrepancy 
with the microscopic counting which gives the finite entropy was resolved by the discovery that the area of the horizon 
is stretched to finite radius once curvature corrections are included. Indeed such corrections have long been known 
to exist in the low-energy effective theories of superstrings [1, 0, [H, [ll], [ll| . The classically computed entropy then 
differs from the Bekenstein-Hawking value lip. UA. lip. Ilal . but agrees (at least up to a coefficient) with microscopic 
countingfpjll [ll M, M, HI HI El H, M, M, fM lli including non-BPS cases [H HH [H, [H, H HE H 13, M, 

Within the models in which the supersymmetric versions of the curvature square terms 
are available, the correspondence was checked using the exact classical solutions [H, [H, [29[. It was also observed 
that good agreement is achieved if the curvature corrections are taken in the form of the Gauss-Bonnet (GB) term 
both in 4D and higher dimensions [13, H [H El • 

To compute the entropy of extremal black holes with the horizon AdS2 x S D ~ 2 from the classical side it is enough 
to construct local solutions in the vicinity of the horizon which is easily done analytically [13, [H, [4{|. But this 
does not guarantee the existence of global asymptotically flat solutions. Construction of solutions with curvature 
corrections, apart from purely perturbative probes I50l. l5ll. 1521. req uires numerical inte grat ion of the field equations. 



M, HI, [M 11, HI llll. The global existence of 



For non-extremal black holes this was done in [53l . [54j, l55l . 
extremal black holes in the 4D model with the GB terms endowed with an arbitrary dilaton coupling a was proven 
in [gl, [65[ . It turned out that global asymptotically flat black holes with the horizon AdS2 x S 2 existed for the dilaton 
coupling below the critical value of the order a CI ~ 1/2 and less than i. This range does not include the heterotic 
string value a = 1 nor i. This result is modified in the presence of the magnetic charge [66j . which extends the 
region of the allowed couplings and serve as the order parameter ensuring continuous transition to the theory without 
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curvature corrections. It is worth noting that our model has neither continuous nor discrete S-duality, so properties 
of the purely electric solution essentially differs from that of dyons. 

The purpose of the present paper is to investigate existence of global solutions for small stretched purely electric black 
holes in higher dimensional Einstein-Maxwell-Dilaton theory with the Gauss-Bonnet term (EMDGB). We construct 
the local solutions in terms of series expansion around the degenerate event horizon for an arbitrary space-time 
dimension D and calculate the discrete sequence of black hole entropies using Sen's entropy function approach. The 
entropy is found to be monotonically increasing with D. Then we continue numerically these local solutions and 
show that in dimensions higher than four the heterotic string value of the dilaton coupling lies inside the range of the 
existence of global asymptotically flat static black holes. 

We also investigate physical significance of the so-called turning points which were encountered in numerical so- 
lutions within the four-dimensional EMDGB theory [53|, H3, [5f| [56|, [64], [6(|. They correspond to mild singularities 
at finite radii outside the horizon where metric and its first derivatives are finite, but the second derivatives diverge. 
Numerical solutions can be extended through these singularities, which we call 'cusps' in this paper, by suitable re- 
definition of the integration variable [56l |69| . In four dimensions the solution extended this way then meets a stronger 
singularity at finite distance, so actually the cusp is just a precursor of the strong singularity. In higher dimensions 
(D > 7) we encounter an interesting new feature: the cusps come out in pairs of right and left turning points, so the 
extended solution finally may be even asymptotically flat. This could correspond to a novel type of black hole coated 
by cusp pairs. But somewhat disappointingly, our analysis shows that continuation of geodesies through the cusp 
singularities in the extended manifolds cannot be performed in a smooth way. Thus we are inclined to reject such 
extended manifolds as physical black hole solutions. Instead, we interpret the occurrence of cusp singularity as failure 
to produce asymptotically flat black holes. This gives an upper bound on the dilaton coupling. We find numerically 
the sequence of critical dilaton couplings for 4 < D < 10 which turns out to be increasing with D. 

Another novel feature of EMDGB black holes with a degenerate horizon in higher dimensions is that the role of 
the GB term in the near critical solutions may still be significant. In four dimensions, as was shown in [64], the 
near-critical solutions saturate the BPS bounds of the corresponding theory without curvature corrections. This 
means that relative contribution of the GB term becomes negligible when the dilaton coupling approaches its upper 
boundary. We find that for D > 7 this is not so, and the BPS conditions are not satisfied in this limit. 

This paper is organized as follows. In Sec. II, we define the action, present the field equations in various forms 
and discuss symmetries of the system. In Sec. Ill, we review solutions for small black holes without GB corrections 
as well and the solutions with the GB term but without dilaton. Then we construct the local series solutions near 
the horizon and calculate the entropy of stretched black holes using Sen's entropy function. We obtain the discrete 
sequence of the entropies of curvature corrected black holes in various dimensions interpolating starting with twice 
the Hawking-Bekenstein value A/2 for D — 4 up to 41A/52 for D = 10. In Sec. IV, we present asymptotic expansions 
of the desired solutions, introduce global charges and discuss the BPS conditions. The next Sec. V is devoted to the 
cusp problem. We explain why extension of solutions through the cusp singularity is physically unacceptable. Finally 
in Sec. VI we present numerical results for various dimensions and explore the fulfillment of the BPS conditions on 
the boundary of the allowed dilaton couplings. 



II. SETUP 

A low-energy bosonic effective action for the heterotic string theory with the curvature corrections is given by [ill . [l2j 

1 = ltc J d ° x ^- Z ^ + *~ 2 ^ M * - + j^GB^j , (1) 

where F^ v is the Maxwell field (we use a truncation involving only one Z7(l) field), £gb is the Euler density 

Cqb = R 2 — 4:R fJiV R ,J "' / + R a f)nvR al3 ^ v ', (2) 

and a' is the Regge slope parameter. The tilde denotes the quantities related to the string frame metric g^ u . The 
action can be transformed to the Einstein frame with metric g M „ by the conformal rescaling 

g^u = g»v, (3) 

giving 



16ttG 



J d D x^-g (r - - Q^F^F^ + ^4>t^£ gb + f(d$, R)^j 



(4) 
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where !F(d$>, R) denotes the cross terms of 9$ and curvature coming from the GB term under the frame transforma- 
tion. For simplicity, we do not include these terms in our analysis. We expect that inclusion of these terms might 
affect the black hole properties only quantitatively but not qualitatively. Then redefining the dilaton field as 

$ = eV 2 (^- 2 ) * (5) 

we obtain the action 

I= T^gJ dDxV ~ 9 i R ~ ~ e 2VJR5 ^ * V" + ^e 2 VV(^ ^£ GB ) . (6) 

In this action we have the sequence of the dilaton couplings 

«L = (7) 

relevant for the string theory. If we do this in 4 dimensions, we have the dilaton coupling a s tr = 1, but if we do this 
in 10 dimensions, we have a s tr = 1/2. It will be convenient, however, to consider the above action for two arbitrary 
dilaton couplings a and b: 

I = I d D x^—g (R - 2d^4> ~ e^F^F^ + ae 2b *£ GB ) , (8) 

where we also denoted the GB coupling a' /8 = a. 

The space-time metric is parametrized by two functions cu(r) and p(r): 

(if 

ds 2 =-u(r)dt 2 + — +p 2 (r)dn 2 D _ 2 . (9) 

For convenience, we list in Appendix A the relevant geometric quantities for more general static spherically symmetric 
metrics. 

We will consider only purely electric static spherically symmetric configurations of the B-dimensional Maxwell field 



Then, integrating the Maxwell equations 



one obtains 



A = -f(r)dt. (10) 



(p D - 2 f'e 2a *y = 0, (11) 



f{r) = q eP 2 - D e- 2a t (12) 



where q e is the electric charge, which is considered as a free parameter (note that the physical electric charge defined 
asymptotically differs from this quantity, see Sec. lIV A|) . 

A. Field equations 

We present the Einstein equations in the form 

G^ = 87rG(T- at +T« B ), (13) 

where T™ at is the matter stress-tensor 



87tGT™ = 2 



(14) 



8^3^ - hcl^^g^ + e 2a * (f^F^ - \F aP F af3 g„ 

and the T^ 8 is the effective gravitational stresses due to the GB term 

8irGT^ B = -ae 2 ^ [H^ + 8 (26 2 V Q 0V /5 (/» + ftVW) P mv p] , (15) 
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where 



H^v — 2{RR pv — 2R fia R a u — 2R ap R liai/ p + R^upyRv" 7 ) — ^gb g^v 

For the metric ([9]), the components of are 

(D - 2)u) 



G u = 

Gaa = 



2p 2 

D-2 



\2ujpp" + puj'p' + (D- 3)(tup' 2 - 1)] 



2up 2 



pjp' + {D-3)(lu P ' 2 -1)] 



* p V + [ 2LU pp" + 2pJp' + (D- 4)(cp' 2 - 1)] , 



while the energy-momentum due to matter fields is given by 

8nGT£ at = ^V 2 + e 2a< W 2 , 

2a<j>J 

) 



8irGT„ 



mat 
rr 



b 12 - e 2a 



8irGT^ at 



-p 2 (u<p 



12 _ 2a</> j/2 



The energy-stress tensor due to the GB term is more complicated 



8ttG 



rpGB _ 



*2b</> tt 



-^W' Wp')(-p' 2 i) - ^W 2 - i) 2 

!^3tL [{2u<t>" + uj'4>'){wp' 2 - 1) + 2up'ct>'(2ijp" + to'p')] 



4bD 2 oj 2 p'(j)' 



(^ ,2 -i), 



8?rG rpGB 

2b<t> rr 



D 2 



a e 



3 v ~'* " ' 2wp 4 ' 



+26 



-( W p' 2 -l) + -^ 7 (u;p' 2 -l) 



-(3wp' 2 - 1) + 



4 W, a 
[up 



1 



8ttG 



rpGB _ 
-L flfl — 



wp p° 

9 D 3 D 3 
D|[ w >p' 2 - 1) + 2^'p'p" + u l2 p' 2 ] + t^L(up')'(Lop' 2 - 1) + ^|(^p' 2 - l) 2 

p 2p 2 



+46 



DIp(ujlj'p'4>')' + Dt{w(t)')'(ujp' 2 - 1) + 2DlLjp'4>'((jp')' 



D 3 5 up'cj>' L _ , 2 



( W p' 2 - 1) 



+ 86 2 [DfWp/oV 2 + L>|w<// 2 (wp' 2 - 1)] , 



where we have introduced the dimension-dependent coefficients 

D% = (D- m)„ = (I? - m)(D - m - 1) • •• (D - n), n > m. 
The dilaton equation reads 



2(w</>')' + 2Dlu<f>'— + 2a/' 2 e 2 
P 



vbD 2 e 2b(1, < 2 



W 2 -l)] ; , n4 (V 2 -!) 2 

~2 h ^5 Ti 



/0 



4£*( W p') 



(16) 
(17) 



(18) 



(19) 



(20) 



(21) 



= 0. (22) 



From the Einstein equation, one can derive the following two second order equations for the metric functions p(r) 
and u>(r), 1 which are more convenient for numerical integration: 



D 2 2 pp" + 2p 2 <// 2 - 4a6£> 2 [(up 12 - l^'e 26 *]' + 2a£ 2 e 2fc * 2bJp l2 $ - D 



1 



P" = o, 



(23) 



1 Namely, the equation for p is — ^- [(Einstein equation)^ + oj 2 (Einstein equation),.,-] and the equation for u) is ^(Einstein equation)^. 
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pw" + 2Dl(ujp')' + Dl Wp,2 ~ 1 + 2ujp<j)' 2 - 2pf' 2 e 2a * - 8abD 3 (uouo' p> $ e 26 *)' 
-aDle 2b ^Dl {uJp ' 2 ~ 1)2 + ADl [(cop')' + 2bup'4>>] 

+2 [u" + 46(w0')' + 8&W 2 ] ~ 1 + 2^ [Wp" + w' 2 // + 8&u^'(o;p')'] [ = °- ( 24 ) 



B. Symmetries of the reduced action 

One can check that equations of motion are invariant under a three-parametric group of global transformations 
which consist of the transformations of the field functions: 

^< P^pe S , 0^0+-, f^fe%-% 5 , (25) 
accompanied by the shift and rescaling of the radial variable 

r -> rc^ +s + v. (26) 
Transformation of the electric potential is equivalent to rescaling of the electric charge 

qe ^ qe e( D - 3+ i) s . (27) 

Not all of these symmetries are the symmetries of the Lagrangian, however. Integrating the action {8} over the 
(D — 2)-dimensional sphere and dropping integration over time integral, one obtains the one-dimensional reduced 
Lagrangian from the relation I = J Ldr. Up to the total derivative one has: 

L = Dip' (up - 3 )' + Dip -* - 2p D - 2 {^ 12 - f' 2 e 2a *) 

- Ud 2 4 p' 3 (co 2 p D - 5 e 2 ^)' + AaD 2 iP ' (^"V^) ' 



[AbD 2 p D - 4 Lo'4>' - 2D 2 p d - 5 lu' P ' - D 2 p D - 6 (cop' 2 - 1)] {up 12 - 1). (28) 
It is easy to check that the one-dimensional action remains invariant under the above transformations provided 

p, = —2(D — 3)5, (29) 
namely under the following two-parametric group of global transformations: 

r^rc-^l,, ^ WC - 2 ( M ) S , P^pc 5 , ^-0+?, f->fe-( D - 3+ *) s . (30) 

o 

They generate two conserved Nocthcr currents 

J ' i= {l&** A - L ) d ' r 



d r Jg = 0, (31) 

3=0 



where Q A stands for w, p, <f>, f, and g = 5, v. The conserved quantity corresponding to v is the Hamiltonian 

H = D 2 2 p' (lo P d - 3 )' - D 2 3 p°- 4 -2u>p D -y 2 + 2p D - 2 /' 2 e 2 ^ 

- ac 2 ^ [4W? 2 p D - V0'(3cV 2 - 1) - 2D 2 p D - 5 u;'p'(3u J p' 2 ~ 1) - D 2 5 p°- 6 (ojp' 2 - l){3ujp' 2 + 1)] 

- AaD 2 p' 3 (u; 2 p D ~ 5 e 2 ^y + 4aD 2 p' (wp^e 2 ^)' . (32) 

This is known to vanish on shell for diffcomorphism invariant theories, H = 0. The Noether current corresponding to 
the parameter 5 leads to the conservation equation d r Js = 0, where 

J 5 = -DjrH - D\J 'p°- 2 + ^lu P d - 2 ^' + 4(D - 3 + q e f 

+ ae 2b * [(cop' 2 - l){2D 2 2 D 2 3 u;'p D - 4 - 8bD 2 3 iop D -^') + 8bD 2 3 uLo' p°- 3 p' cj)'} . (33) 

Symmetry transformations will be used to rescale numerically obtained solutions to desired asymptotic form and 
obtain true physical parameters of the solution. 
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III. STRETCHING THE HORIZON OF SMALL BLACK HOLE 
A. Small dilatonic D-dimensional black hole without GB term 

Let us first discuss the black hole solution without the GB term. It can be presented in the form (6?| 

ds 2 = -f + fl 1+ $^dt 2 + f+ 1 f_ l+1 ^~ TU ^dr 2 + r 2 f^~ T ™dn 2 D _ 2 , (34) 

c 2a * = e 2a ^ fZ^ , F tr = UD - 3) [r+r -J_ 2 e -^°°^^, (35) 

VA r u z 

where 



The mass and the electric and dilaton charges are given by 



M = 



16ttG 



(D - 2)(r^ 3 - r D - 3 ) + 4 ^ A 3) r°' 3 



(37) 



„ _ (39) 

For a = 0, this solution reduces to the D-dimensional Reissner-Nordstrom solution. In the extremal limit r + = 
r_ = ro, it contracts to 

ds 2 = ~f 2 dt 2 + f^dr 2 + r 2 dn 2 D _ 2 , / = 1 - ^J, (40) 

and has a degenerate event horizon Arf52 x S D ~ 2 . Note that for a = 0, the GB term decouples from the system, so 
this solution remains valid in the full theory with a^O. 
For a ^ 0, the extremal solution reads 

4QD-3) 2(J>-4) 4 2 4 

^ 2 =-/o (£ " 2,A * 2 + /o °~ 3 (D - 2)& dr 2 + r 2 f D - 3 (D - 2)& dtf D _ 2 . (41) 
This has a null singularity at the horizon. The Ricci scalar in the vicinity of this point diverges as 

R~ ( r ^-3 _ r( D-3)-u 2 r3 + (D _ 4 2)Aj (42) 

together with the dilaton function 

e 2 O0 ^ (r o-3 _ r ^-3)-¥. (43) 
The divergence of the GB term near the horizon is 2 

e 2a0 £ G s|r=r + - (r+ - r-)"^^-^^), (44) 

so one can expect that the GB term will substantially modify the dilaton black hole solution in the extremal limit. 
The mass, the dilaton charge and the electric charge for this solution (defined as in Sec. IV below) are 

_ 47rG A ' ye _ Qe * ~ 4ttG VA ' _ 47rG A (45j 

They are determined by a single parameter ro, so we have the following relations among the three quantities 

V = aM, Q e = y/AM, (46) 

which imply the following BPS condition 

a 2 M 2 + V 2 = ^-Q 2 e . (47) 



2 We use this occasion to correct Eq. (33) of our previous paper for D = 4 64j. 
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B. Wiltshire black hole 

Another limit in which our action admits an exact solution is that of vanishing dilaton. This is consistent with the 
field equations for a = b = 0. In this case an exact solution was found by Wiltshire [HI] 



, , 1 r 2 I I UirDlaM 8Djaq 2 \ 
^ = 1+ 2Dfr[ 1 *f + W a r°-x " DlA^V ) ' (48) 

The lower sign corresponds to an asymptotically AdS space-time for a > and to an asymptotically de Sitter solution 
for a < 0. The upper sign leads to an asymptotically flat solution coinciding with the D-dimensional Reissner- 
Nordstrom solution. These solutions exist in dimensions D > 5 where the GB term is not the total derivative. The 
asymptotically flat solution has two horizons which coincide in the extremal limit for a special value of the electric 
charge. For the extremal solution, the mass and the charge can be expressed in terms of the single parameter, the 
radius of the horizon tq: 

M = ^f(D - 2)[rg + (D- Afa]r%-\ <ll = ^H + Dfa]^. (49) 
Conversely the radius can be expressed as 




An DIM I ( A-itD\M \ 2 2Djq 2 



£>|o D _ 2 1 V \D 2 n D _ 2 ) d 2 



(50) 



C. Local solution near the horizon 



In what follows we set b = a as relevant for the heterotic string theory case, but still keeping a arbitrary. Assuming 
that the full system with the GB term admits the AdS% x S D ~ 2 horizon, r = rjj , we look for the series expansions of 
the metric function in powers of x = r — rn : 



PiX 



P(r) := e 



2a4>(r) 



i=2 



i=0 



5> 

i=0 



(51) 



The function ui starts with the quadratic term in view of the degeneracy of the horizon, while two other functions have 
the general Taylor's expansions. Denoting the physical radius of the horizon po = p(rn), we obtain for the leading 
order coefficients: 



U>2 = 



Po = 



Pi 



Aa(2D - 7) : 



and po is related to the electric charge via 



4VM2D - 7) 
y/Dl{D*-D-i 



(52) 



(53) 



Note that the expression under the square root and the right hand as a whole are positive for D > 4. The horizon 
radius is fixed entirely by the electric charge, like in the extremal Reissner-Nordstrom case. In our units the GB 
parameter a has dimension L 2 . When the GB term is switched off (a — > 0), the horizon radius shrinks, as expected 
for small extremal black holes. Higher order expansion coefficients exhibit dependence on only one free parameter, 
namely the P\ in the dilaton expansion. Other coefficients are expressed in terms of the horizon radius po and P±, 
the first sub-leading coefficients being 



u 3 



Pi = 



2<xPi 



Dl(3D 5 - 43L> 4 + 213L> 3 - 421D 2 + 236D + 76)a 2 



3a 2 p^(D 3 -9D 2 + 16D + 

+D\{2D - 7) (3D 4 - 25D 3 + 70D 2 - 56D - 40) 

4<xPi[(L> 4 - 13D 3 + 54D 2 - 72D - 2)a 2 + (2D - 7)(D 2 - 3D - 2)] 
a 2 p (D 3 -9D 2 + 16D + 8) ' 



(54) 
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One can notice that the free parameter enters the expansion coefficients always in the combination Pi /a 2 . This 
facilitates transition to the Wiltshire case. 

In the limit of decoupled dilaton a = 0, the parameter Pj — > 0, while the ratio Pi /a 2 remains finite. In this case 
we have nonvanishing coefficients Po, pi and uf. 

P(r)=P , p(r)= Po +pi{r-r ). (55) 

The asymptotic flatness requires p\ = 1, so we have po = tq. For the Wiltshire solution Po = 1, and we obtain the 
following relation in the extremal case: 

p 2 Q = r 2 =4a(2D-7), (56) 

Substituting KSBJ) into (gHD , we find 




and so our solution with the decoupled dilaton coincides with the extremal case of the Wiltshire solution. 

We can consider subgroup of global symmetry transformation defined by two parameters 6 and p. We can eliminate 
the parameter po from the expansion on the horizon if we apply the transformation with parameters p = — 21npo 
and 6 = — lnpo- The other transformation with parameters p — 21n|Pi|, S = can take out Pi from expansions. 
Choosing the absolute value |Pi| in the second transformation allows us to get the remaining parameter £ = jp^j in the 
expansion at the horizon which fixes the sign of p (pi > to obtain a global solution) in the expansion. As a result, 
we have a map between parameters of expansion Pi, po and parameters of global transformation p, 6. Typically one 
can first investigate special solution with a simple choice of near horizon data, such as Pi = 1 and po = 1, which are 
the values we use for our numerical analysis below. Then general solutions with arbitrary values of free parameters 
can be simply obtained by the global transformation with p — 2 In 8 = — lnpo- Also it is clear from the relation 
between p and q e that electrical charge plays the role of rescaling parameter. Finally the free parameter Pi could 
be fixed in accordance with the boundary condition at infinity. We can also eliminate the GB coupling constant a 
from the system by introducing new dilaton function F — a P and rescaling charge q e . In this way, we have only two 
parameters, the number of the dimension D and the dilaton coupling a, which affect the dynamics of solutions. 

The values of the integrals of motion (|32| and ([33| in terms of the parameters of the local solution are 

H = -ap°- 2 D 2 P -D 2 p°- 4 + ^p 2 - D , J s = 4(D-2)qJ . (58) 



D. The entropy 

Knowledge of the local solution near the horizon is enough to calculate the entropy of the black hole, assuming that 
the local solution can be extended to infinity. To compute the entropy, we apply Sen's entropy function approach Q 
which is valid for the black holes with near horizon geometry of AdS 2 x S D ~ 2 . Using the notation of [3|] we parametrize 
the near horizon geometry by two constants, V\ and vi related to the radii of AdS^ and S D ~ 2 , as 

ds 2 = vi (-r 2 dT 2 + + v 2 dtt 2 D _ 2 . (59) 



The scalar curvature and the GB term will read 

2D 2 „ D 2 AD 2 , 
R= + -^, C GB = ^t -■ 60 

Vi V 2 U| ViV 2 

The dilaton field and gauge field strength are constant on the horizon 

tf> = u, F Tr = p- (61) 

Sen's entropy function is defined to be the integrand of the action after integrating all angular coordinates of S D ~ 2 . 
Using we obtain 



9 



The parameters vi,V2,u,e are related to the near horizon expansion coefficients by (note the rescaling of time coor- 
dinates T = LL>2t) 

<^ = ^, Po = V^, Po = e 2< ™, q e =pv- 1 v 2 2 ^ e 2au . (63) 

According to the equations of motion, the value of parameters should minimize the entropy function: 

d v J = 0, d v J = 0, d u f = 0, (64) 

which lead to the following constraints 

D 2 A(2D - 7)p 2 2au 2 2(D 2 - D-8) 
V2 = — Vl , v 1= D2 _ D 8 e , P = a, (65) 

and furthermore imply / = 0. These three constraints are exactly identical with the relations (|B"!?|) and (fS"3")) from 
the near horizon analysis. The physical electric charge, q (i.e. Q e defined in subsection IIV A| ) . can be obtained via 
Q = d e f 

q = ^G pV ^ V - 6 = ~teG q - (66) 
The entropy of black holes is related to the entropy function by a Legendre transformation 

S = 2n(qp -f) = 2nqp = S^j^J • (67) 

The horizon area of AdS2 x S 2 is A — vol(f2£)_2)i> 2 2 , thus the entropy can be expressed in terms of area of horizon 

as 

D 3 -D-8. A D 2 -5D + 6 
S - 8(2D - 7)G A -AG + 8{2D -7)G^ SbH + 5gB ' (68) 

and the deviation of the entropy from Bekenstein-Hawking relation by the GB term increases for higher and higher 
dimensions. For example, the ratio of Sgb/Sbh from D = 4 to 10 is 

S GB J 6 10 5 21 281 

^ = l 1 ' 1 '5'y'3'TT'T3/- (69) 

A g eneral discussion on the entropy of theories with quadratic curvature correction and Lovelock theory is given in 



IV. ASYMPTOTICS 



Now consider the asymptotic expansions of the metric function by substituting the following expansions into the 
equations of motion: 



«(r) = l + E§' PM=r + £§. *M = *»+£$ 



(70) 



»=i 



According to the falloff of the Newton potential in different dimensions, one has the first non-zero term in the expansion 
for w and that for dilaton starting from i = D — 3, while p differs from r in (2D — 7)-th terms: 



1 + 



ujd-3 2q 2 e 



-2a0 o 



l 



„D-3 



D 2 r 2(D-3) 



o 



1 



„2D-4 



M P - 3)<? D - 3 1 

p(r) = r — 



(r) = 4>oo + 



(U-2)(2D-7) r 2C - 7 

D-3 1 



»2Z>- 



(71) 



„D-3 



fl? 2 e -2a0 : 

(£>-3) s 



+ ^D-3(PD-3 



1 

r 2(D-3) 



„2D-4 
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One can notice, that these terms of expansion do not contain the GB coupling a. The contribution of the GB term 
is manifest in the third non-vanishing coefficient in p. If the GB term is switched off a = 0, the third non-vanishing 
coefficient of p is 



~ P3D - W 3Dl(3D-10f D - 3 



(D - 3) + aq 2 e e~ 2Q M . (72) 



In presence of the GB term it is 



P2Z3-5 = 2 ^ D -l »euD D - 3 ^ D - 3 e 2a ^ . (73) 



So in D = 4 the GR contribution dominates appearing as p 2 , and in D = 5 both GR and GB contributions appear in 
p~5, but in higher dimensions, GB contribution is leading. 

For the asymptotically flat geometry the global physical quantities, such as mass and charges, can be read out from 
the asymptotic expansion. Since the first sub-leading coefficients are independent of the GB coupling, we can still use 
the formula of global charges for the theories without higher curvature corrections. 



A. Global charges 

The ADM mass is given in our notation by 3 



r B-3( _L _P) _ r D-2 



(?)' 



(74) 



and reduces to 



J M = -^g(£>-2)[w D _ 3 -2(D-4)p 2 ,_ 4 ]. (75) 

For D > 4, in general, the ADM mass could depend not only on the first sub-leading coefficient ujd-3 of ui, but also 
on the sub-leading coefficient pd-4 of p. But we have seen this coefficient is zero, so we have 

WO - 3 = -(D-2)fi D _ 2 - (?6) 

The definition of the dilaton charge V is 



which has a contribution from the expansion coefficient 

A-kGV 



>D-3 



(£>-3)Qz3_ 2 

The physical electric charge can be computed by the flux 



(78) 



Qe = — ^ / cM D _ 2 r D ~ 2 e 2a *~ F tr , (79) 

47TG Jr^oc 

so we have the following relation between this quantity and the charge introduced as an integration constant in the 
previous section: 

Qe = 75— "Qe- (80) 
"D-2 



3 The volume of S°- 2 is Q D _ 2 = ffu^ry and thc 

gamma function is cither T(n + 1) = n! or r(f + 1) = V^Stt for integer n. This 
gives Q D „ 2 = {47r,27r 2 ,|7r 2 ,7r 3 ,lf7r3 2 ,i7r 4 ,^7r 4 } for D = 4, ■ ■ ■ ,10. 
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The asymptotic values of two integrals of motion and ([55]) are 

= [Dl{rp')^ - aDlP^rp')^- 6 ] (uj^p'J - 1), (81) 

Poo 

+4£>|ee/ao - 25^^ P r 2 - 4^ Woo ^-2^. (82 ) 

From H = 0, we can see that Woop'^, — * 1 for r — ► oo, which also regularizes the second integral of motion. For 
Minkowski space (cJoo = p'^ = 1), we have 

=4D 2 qJ oa -2D 2 D 2 M-4Dl-. (83) 

a 

It is possible to apply global transformation to satisfy asymptotically flat condition which fixes one of the free 
parameters in expansion around horizon, P± . Note that the values of the integral of motion are four times of what we 
have in [64j : 

H = \ {uoop'l - 1) , J S = 2q e f oc -M--. 
2 a 



B. BPS condition 



The theory we are considering here does not necessarily have an underlying supersymmetry. However, it is instruc- 
tive to investigate the fulfilment of the no-force condition which is usually associated with the supersymmetry. In 
particular, in the D = 4 case the supersummetric embedding into the heterotic string theory in the supergravity limit 
gives the BPS condition for the extremal small black holes (Q e = Q e e~ a ^°°) 

M 2 +V 2 = Q 2 e . (84) 

This corresponds to vanishing of the sum of the gravitational and dilaton attractive forces and the electric repulsion 



This does not hold if the GB term is turned on. However, it was demonstrated in [64| that on the boundary of the 
allowed domain of the dilaton coupling the role of the GB term is diminished, and the BPS condition is restored. Our 
aim here is to confirm this property. 

In the higher-dimensional cases the gravitational, Coulomb and dilaton forces are 



F 
F 



8irG{D-3) M" 1 
" ^ ~(£>-2)fi D _ 2 rD- 
AttG Ql 



AttG V 2 

n^~2^ 2 ' 



AttG V 2 , . 

F * ~ n=9> ( 85 ) 



so the no force condition reads 



2(D - 3)M 2 + (D - 2)V 2 = (D~2)Q 2 e . (86) 



In the case that the GB term is decoupled, i.e. a — 0, the no-force condition (|86|l at infinity is equivalent to the 
degenerated horizon obtained for the exact extremal dilatonic black hole solutions in the previous section 

V = aM, Q e = VAM =>■ a 2 M 2 + V 2 = — Q 2 e . (87) 

Note that the relation (|86"|) does not involve explicitly the dilaton coupling (though it appears in the definition of 
Q e ). The special case of D = 4 (A = a 2 + 1) was earlier discussed in [6J| in which case <I>i — —2M., 4>i = —T> (with 
different sign convention). 
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V. CUSPS 



It was discovered in [64( that the 4D EGBD static spherically symmetric gravity typically develops cusps at some 
points r = r c in the vicinity of where the metric functions vanish. There they have Taylor expansions in terms of 

V=\r-r c \. (88) 

The metric and its first derivative are regular there, while the second derivatives diverge as y~ x ^ 2 ■ There are therefore 
the cusp hypersurfaces which are the spheres S D ~ 2 of finite radius. These cusp spheres have curvature singularity 
which is rather mild (the Ricci scalar diverges only as y~ 1 ^ 2 , and the Kretchmann scalar as 1/y). They are in fact 
the singular turning points of the radial variable p(r). 

The presence of turning points in the numerical solutions was encountered in the case D = 4 in [55|, HH, H^] . The 
numerical solution can be extended through these points using the technique of [69j and then the solution evolves into 
a strong singularity. Here we find that the situation is similar in higher dimensions D < 6, but starting from D = 7 
the solution can be extended to an asymptotically flat one. 



A. Expansion near the turning points 

The general property of the turning points is that the metric functions and the exponential of the dilaton field, 
/ = {ui(r), p(r), F(r) — ae 2a ^ r '}, have finite first derivative and divergent second derivative, i.e. 

f'(r tp ) = constant, /"(r tp ) -> oo. (89) 

The metric functions and the dilaton can be expanded in terms of the fractional powers of the variable y 

/(«) = /<> ( 9 °) 

i=2 

where we have either y — r tp — r (the right turning point) or y = r — r tp (the left turning point). These two types of 
turning points have opposite signs of the odd-order derivatives, 

/(2«+i) (r _ rtp ) = _ / (2n+l) (rtp _ r)) (91) 

and the expansion coefficients, {fi}, have the same "iterative" relations for both type turning points. The expansions 
read 

uj = Lu Q +aj 2 y + ^3y^ +0(y 2 ), (92) 

P = Pa +P2y + p 3 y^ +0(y 2 ), (93) 

F = F Q + F 2 y + F 3 y? +0{y 2 ). (94) 

They contain four free parameters, namely uq, po, Fq and p 2 (for the fixed charge parameter q e ), other coefficients 
depending on them. The coefficient p 3 is given by the square roots of a second order equation which can have two 
branches (positive and negative) corresponding to double valued solution near turning points. Similarly, and F$ 
also have two-branch solutions. 

The exponents in the turning point expansions are independent of D. Therefore, the rate of divergence of geometric 
quantities is universal. More precisely, the scalar curvature is 

fl „ _3q, 3 po + 2(fl- 2)0^3 1/a 
4 pa 

and the matter stress tensor is finite. Indeed, uj 3 is proportional to p 3 which has double values (with opposite sign) 
near the turning point. Therefore, the sign of divergent scalar curvature also changes. Moreover, one expects that 
the GB combination should have y~ x divergence, but actually it is weaker, namely y~ x / 2 . 

For numerical integration, we rewrite the equations of motion as a matrix equation of the dynamical system 



Ax' = b, 



(96) 
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where 6D vector x denotes x(r) — {uj(r), uj'(r), p(r), p' (r), F(r), F'(r)}. The solution is ill-defined at the points where 
det A = 0. The turning points are special cases of the general situation (see [69] for complete classification). We can 
extend solutions through the turning points introducing a suitable new parameter a: 

(IT 

r= — = A det A, (97) 
da 

and generalizing the dynamical system to one dimension more, i.e. x(a) = {uj(a),uj(a), p(a), p{a), F(a), F(a),r(a)}. 
The matrix equation then becomes 

A * = h > ^=(o"l 6 )< ^{xIa)' ( 98 ) 

The parameter A can be fixed by normalization of x. We choose x 2 = x 2 + f 2 = 1 and this ensures that r is finite for 
both small and large values of det A which is useful for numerical calculation. 
In terms of a, we have the following result near turning point 

r{o) = rtp + n (o-tp - ct) 2 + • • • , (99) 

where r\ < for the right turning points and r\ > for the left ones. The metric can be rewritten as 

ds 2 = -Lj(a)dt 2 + + p 2 (a)dtl 2 D _ 2 , (100) 

where W = u>/r 2 . Now, the functions u>, W, p are single valued functions of a. 



B. Geodesies near the turning points 

As mentioned before, the curvature weakly diverges at turning point. One would expect the property near turning 
points is much better than near the singularity. So let us check the radial geodesic of t and a as functions of the 
proper time A (d9/dX = = d<p/dX). The relevant Christoffel symbols are 



1 ta — „ j 

2 uj 



Iwuj, 



1W 

2W 7 



(101) 



The geodesic equation for t 



dH 
d\ 2 



uj dt da 
uj dX dX 



(102) 



can be simplified as 



d_ 

dX 



dt 
dX 



(103) 



or 



dt 
dX 



C 

* 



(104) 



where the integration constant C > means the "energy" per unit mass of test particle at infinity. The geodesic 
equation for a coordinate is 



fa 1 

dX 2 + 2 



dty w /day 
Wuj[ dx) ~w{dx) 



= 0. 



(105) 



After integration, it reduces to 
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or 



fdaY (C 2 \ C 2 -kuj 



= W\—-k \ = —72—, (107) 



where k = 1 for time-like geodesic and k — for null geodesic. 
The geodesic solutions are 

t=-A, (a-a tp r = ~^X 2 ± VC ]- kLV0 X. (108) 

There are two possible solutions for ex, the minus branch is valid for — oo < A < and plus branch for < A < oo and 
both geodesies are terminated at the turning point (A = 0). The only possible extension for the geodesic solution is 
"gluing" these two solutions, i.e. 



(CT-ctp) = -^a X +sign(A) — A. (109) 



However, one can easy check that the result is not a smooth solution at A = and the second derivative of a with 
respect to A generates a delta function. Therefore, such extension, in general, is not a solution of (|105[) . Similar 
situation happens for the time-like geodesic. Hence, the cusp turning points are not extendable for the geodesies. 
However there is an exception for the special values of C 2 = fcc^o and ku>2 < (if u>2 > the time-like geodesic cannot 
reach the turning point). 

It is instructive to compare divergences in various cases. To study this, we use Krechman scalar K = R a ^ s R a ^js- 
The Schwarzschild black hole at r = has K ~ r~ 2 ( D ~ x > , anc j the Reissner-Nordstrom black hole at r = has 
K ~ r ~ 4 ( D ~ 2 \ GB pure gravity black hole (D > 5) at r — has 1) K ~ r -2 ' 23-1 ) for charged solution, and 2) neutral 
one K ~ r^- 0-1 ). The extremal dilatonic black hole with GB term at a turning point K ~ for any dimension: 

R = ± ^2 WoP3 (D-2) +PoW ^ 

where y — r tp — r for upper sign (a first type of turning point) and y = r — r tp for lower (a second type). If p 3 = 0, 
the scalar curvature at such a type of point is regular, but it imposes some constraint on the parameters because p 3 
depends on the parameters. 



VI. NUMERICAL RESULTS 

Lets consider the special limit a — > in which the dilaton field decouples. The analytical general solution is 



"V - 1+ 2l^{^f + (D- 2)v l(^ 2 )^ - 15pA + ) - ^ 

For the case A = 0, the radius of degenerated (extremal) horizon is r 2 , — A(2D — 7)a which can be obtained from (|56p . 
It is clear that the horizon shrinks to a point when we turn off the GB term. Although the dilaton field is decoupled, 
in the dimensions D > 5, the GB term still gives a non-trivial contribution to equations of motion which will break 
the BPS (non-force) condition. In more detail, the ratio of the mass and the electrical charge for an extremal solution 
can be computed 

^L = ^ 12 ) 2 >! (112 ) 

Ql 4(D 2 -D-8)(2D-7) " ' ( ' 

and our numerical analysis give consistent results in the decoupling limit. Moreover, the numerical analysis also 
indicates that the dilaton charge is proportional to the dilation coupling times mass and 

a 2 M 2 , , 

-p5- ^ L ( 113 ) 

but the exact form for the ratio is still unknown. In these two ratios, the equality holds for D = 4 which saturates 
the BSP condition ([57]). 
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In the numerical results, we are going to show the following quantities. From the symmetries (|30l) . we know that if 
we regard the electrical charge as a scaling parameter, the following ratios will depend only on the dilaton coupling: 



kM(a) = — , k D (a) = — — , k F (a) = . (114) 

v^e t^:e 



For verifying the BPS conditions, we will analyze the ratios 



a 2 M 2 AM 2 , A(a 2 M 2 +V 2 ) . , 

k BPS = TTTTr? ■ V 115 ) 



V 2 ' Ql ' 2a 2 Q 

We now present our numerical results. 



A. D = 4 



For convenience, we first recall the results for D — 4 found in (64|. Starting with small a (for a — one has the 
Reissner-Nordstrom solution) one finds that asymptotically flat black holes with degenerate event horizons exist up 
to a — a CI = 0.488219703. The critical value of dilaton coupling a cr separates the regions where there exist regular 
asymptotically flat solutions for a < a cr and the singular ones for a > a cr , which firstly have a throat (p 1 = 0), then 
a turning point where p" changes sign and p" diverges and finally a singular point. In the limit a — > a cr , the mass 
A4 diverges, and somewhat surprisingly, the BPS condition of the theory without curvature corrections holds for the 
ratios of parameters. This can be understood as dominance of the Einstein term over Gauss-Bonnet. Indeed, if we 
keep the mass fixed, the limit corresponds to gravitational constant G going to zero. Then the Einstein term becomes 
greater, unless the GB term is increasing similarly. 

The critical value of the dilaton coupling is less than the heterotic string value a = 1 or 1/2, so no asymptotically 
flat extremal EMDGB black holes exist in D = 4. 



B. D = 5 and 6 




FIG. 1: Left: fe_M(a),fco(a) and kp(a) (multiplied by a factor 10 3 ) as functions of a in D = 5. Right: ksps and the ratios 

AM 2 a 2 JVl 2 



The critical value of dilaton coupling for D — 5 is found to be a^ -5 = 0.872422988. The corresponding string value 
([7]) is smaller: a^, =5 = 0.816496580, so asymptotically flat stretched dilatonic black holes exist in five dimensions. 

For D = 6 we obtain a°= 6 = 1.432881972, the corresponding string value Q being a°= 6 = 0.707106781). We 
observe that while the critical dilaton coupling is increasing with D, the string value ([7]) is decreasing, so we expect 
that for D > 4 we will have always a str < a cr ). 
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The metric functions and the dilaton exponential F for asymptotically flat black holes are given in Fig.[5]for a = 0.3 
and 0.8 (which are smaller than the critical value in D = 5), and in Fig. [5]for a — 0.4 and 1.4 in D = 6. Those for a 
larger than the critical value are displayed in Fig. [3] in D = 5 and in Fig. [5] in D = 6. 




FIG. 2: Radial dependence of metric functions ui, p' and the dilaton exponential F for asymptotically flat black holes (a < a£ 5 ) 
with dilaton couplings a\ = 0.3 and ai — 0.8 in D = 5. 




FIG. 3: Radial dependence of metric functions u>,p' for singular solutions (a > a® 5 ) for dilaton couplings a — 0.9 in D — 5. 
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FIG. 5: Radial dependence of metric functions ui,p' and an exponential of dilaton field F for asymptotical flat black holes 
(a < a£ =B ) of dilaton couplings ai = 0.4, ai — 1.4 in D = 6. 




FIG. 6: Radial dependence of metric functions u>,p' for singular solutions (a > a£ B ) of dilaton coupling a = 1.5 in D — 6. 

C. D = 7 

The critical value for the dilaton coupling a^. =7 = 1.793909999, the heterotic string value being a®= 7 = 0.632455532. 
As in lower dimensions, the critical value corresponds to the appearance of the first cusp (turning point) in the solution. 
The novel feature in D = 7 is that after the right turning point the left one appears, and the solution can be extended 
along the lines of 56]. Using the same procedure one can extend the solution to an asymptotically flat one, as 
shown in Fig. \7\ With further increasing dilaton coupling the number of pairs of the turning points increases, so the 
asymptotically flat extended solutions look as shown in Fig. [8] The global parameters change in step-function-likc 
manner each time when one new turning point is created, see Fig. [5] However, the extended solution cannot be 
considered as true black hole solution, since geodesies, as we have shown in Sec. V B, cannot be continued smoothly 
through the cusp singularities. Therefore we have to consider the critical value of the dilaton coupling in D — 7 as 
the true boundary of the range of a. 

In the case D = 4 [64] it was observed that on both boundary of the dilaton coupling a — ► and a — > a cr , the BPS 
condition of the EMD theory is saturated. This can be understood as indication that the GB term is decoupled in 
these two limits. Indeed, for a = it is obvious, while in the limit a — ► a cr the mass tends to infinity, in which case 
the Einstein term turns out to be dominant. In higher dimensions situation is different. For D > 5 decoupling of the 
dilaton a = does not switch off the GB term; instead we have to deal with Wiltshire solutions of the EMGB theory. 
So the BPS saturation is not expected for a = 0, and this is confirmed by numerical calculations. For D — 5,6 the 
BPS condition still holds on the right boundary of the dilaton coupling a — > a CI (see Figs, [l] [4|. However, in D = 7, 
the BPS condition does not hold anymore (see Fig. [TU)) . This means that the GB term does not decouple in the limit 
a — > a cr as in lower dimensions. 
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FIG. 7: Dependence of metric functions w and p' on radial coordinate for a=1.8inD = 7 with two turning points. (Insets 
describe properties of solution nearby turning points). 




FIG. 8: Dependence metric functions w and p' of radial coordinate for a = 1.81 in D — 7. There are fourteen turning points. 

D. £> = 8,9,10 

Properties similar to those in the case D = 7 were observed for higher dimensional solutions D — 8, 9, 10. The critical 
values of dilaton coupling are a£= 8 = 1.887653885 (a^= 8 = 0.577350269), af r = 9 = 2.002906751 (a^= 9 = 0.534522483) 
and a°= 1Q = 2.121748877 (ag= 10 = 0.5). The numerical results are presented in Figs. [H| [TBI TO TO ITS] TO The 
BPS condition is not fulfilled on both boundaries of a. Supercritical solutions can be formally continued to infinity 
(as asymptotically flat) through the cusps which are met in pairs like in the case of D = 7. However we do not qualify 
them as physical black holes for the reasons explained in Sec. V. 

VII. CONCLUSIONS 

Here we summarize our findings. First, we have constructed explicit local solution of the EMDGB static extremal 
black holes in the vicinity of the horizon and calculated the corresponding entropies. The ratios of entropy to the 
Hawking-Bekenstein one A/4 increases from 1 for D = 4, 5 to 41/13 for D = 10. The entropy does not depend on 
the dilaton coupling. Contrary to this, the asymptotic behavior of the solutions crucially depend on dilaton coupling 
and asymptotically flat black holes exist only for a < a cr . The critical value of the dilaton coupling depend on D 
and increases with D. For D = 4, a cr is smaller than the heterotic string value, therefore no stretched black holes 
exist in the effective heterotic theory. In contrast, for D > 5 the heterotic values of a lie inside the allowed region. 
Numerical solutions for asymptotically flat black holes are constructed for 4 < D < 10. We investigated the ratios of 
the mass, dilaton charge and electric charges which show the degree of deviation form the BPS bounds in the absence 
of GB term as functions of the dilaton coupling. It is observed that for D < 5 the BPS bound is saturated near the 
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FIG. 9: Left: kM(a),kD(a) and kp{a) in the region before formation of turning point in D = 7. Right: same quantities after 
formation of turning points (the number of turning points is denoted by italic font). 




FIG. 10: Left: Ratios of ksps and , a in the region before formation of turning point in D = 7. Right: those after 

formation of turning points (the number of turning points is denoted by italic font). 

threshold value a cr , thus demonstrating that the contribution of the GB term is effectively small there. For larger D 
such a behavior was not observed, indicating that the GB term remains important on the boundary. 

The failure to reach the flat asymptotic in numerical integration manifests itself as emergence of turning points of the 
radial variable in which the scalar curvature has very mild divergence. The solutions then exhibit typical cusp-shaped 
behavior. It was suggested before that these turning points should be passed by changing the integration variable in 
a suitable way so that the solution can be continued through these singularities. We have found that in dimensions 
D > 7 the turning points comes in pairs, and the solution can be formally extended to the flat asymptotic. However 
an inspection of radial geodesies reveals that they cannot be analytically continued through cusp singularities, so we 
do not believe that continuation of numerical solutions through the cusps is physically meaningful. 
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FIG. 11: Dependence metric functions w and p' of parameter a for a — 1.8 in D — 7 with two turning points, (insets describe 
properties of solution nearby turning points) . 
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FIG. 12: Dependence metric functions w and p' of parameter a for a — 1.81 in D — 7 with fourteen turning points. 



APPENDIX A: GEOMETRIC QUANTITIES FOR SPHERICAL SYMMETRIC METRIC 

This appendix gives detail geometric quantities associated with the following spherical symmetric metric in and 
dimensions D 

ds 2 = _ e Mr) dt 2 + e Mr) dr 2 + Q 2 w (r) ^ ( A l) 

where k denotes the spatial curvature. The Riemann and Ricci tensors have the following components 

R trtr = e 2u (u" + u' 2 -uV), 
Rtatb = e 2u - 2v+2w u'w' g ab , 

Rrarb = ~C 2w ( W " + w' 2 - v'w') 9abl 

Racbd = (-e 4w - 2v w' 2 + kc 2w )(l-S ac S bd )g ab g cd , 

R tt = c 2 ^(u" + u'H'), 

R rr = -(u" + u' 2 -u'v')-(D-2)(w" + w' 2 -v'w'), 

R ab = [-e 2w - 2v (w" +w'H') + k(D-3)]g ab , (A2) 
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FIG. 14: Left: kM{a),kD(a) and kp(a) in D = 8 in the region before formation of turning point. Right: those after formation 
of turning points. 



and then the scalar curvature, Ricci square {R 2 ^ = R^W) and Riemann square (i?^ = Rafi^vR^^) are 

R = -e' 2v [2u" + u'H' + u' 2 - u'v' + D 2 (2w" + w'H' + w' 2 - v'w')] + kD 2 e~ 2w , 
= -c- 2v [2(u" + u' 2 - u'v') + 2D 2 (w" + u'w' - v'w' + w' 2 ) + D 2 (w' 2 - kc 2v - 2w )] , 
R\ v = c- 4v (u" + u'H') 2 +c- 4 >" + u' 2 -u'v' + D 2 (w" + w' 2 -v'w')} 2 

+ D 2 [c- 2v (w" +w'H')-kD 3 3 c- 2w ] 2 , (A3) 

l a/3ftv 

+ 2D 2 3 (e~ 2v w' 2 -fee- 
where H is defined 

H = u — v + (D — 2)w, (A4) 

and the following notation is used 

D™ = (D-m) n = (D -m)(D-m- 1) • •• (D -n), n > m. (A5) 



K&» v = 4 c ~ 4v («" + u ' 2 - u ' v ' f + 4 ^2 e" 4 W 2 + AD 2 e~ 4v (w" + w' 2 - v'w') 2 
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FIG. 15: Left: Ratios of k BPS and ^-,2^- m D 
formation of turning points. 



8 in the region before formation of turning point. Right: those after 
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FIG. 16: Left: kM(a), fcu(a) and kp{a) in D = 9 in the region before formation of turning point. Right: those after formation 
of turning points. 

The GB combination is 

£ GB = D 2 s e- 4v ^Di( W ' 2 -ke 2v - 2w ) 2 + 8u'w'(w"+w' 2 - W 'v') 

+ 4 [u" + u' 2 - uV + Dj(w" + u'w' - v'w' + w' 2 )] (w' 2 - k e 2 ^ 2 ™) J 



AD 2 e~ u - v - 2w [c u - 3v+2w u' (w' 2 -ke 2v ~ 2w )]' 



+ AD 2 c- iv {w" + u'w' - v'w' + w' 2 ) (w' 2 - k e 2v - 2w ) 
+ D 2 e- 4v (w' 2 -ki 



>r t ..,'2 u e "2v-2w\ 2 



(A6) 



One can easy check that, in four dimension, the GB term of ^/^gCcB is a total derivative. For the gauge choice of 
coordinates 



u = — v = ^ In oj, w = In p, 



(A7) 
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FIG. 17: Left: Ratios of k BPS and in L> = 9 in the region before formation of turning point. Right: those after 

formation of turning points. 
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FIG. 18: Left: kM{a),kn{a.) and fa?(a) in D = 10 in the region before formation of turning point. Right: those after formation 
of turning points. 



the relevant quantities become 

R = -p- 2 [uj" p 2 + 2D 2 2 p{Ljp" +lj'p') + Dl(up' 2 -k)] , 

„2-D 



P 



iU p 



'» D - 2 + 2D 2 ujp'p D - 3 y + D 2 p'(up D - 3 )' + D 2 p D - 4 k 



Cqb - D 2 3 p- 4 {2uj'p'p 2 {1ujp" + u'p') + 2 P [lo"p + 2D\{up')'\(iop 12 - k) + D\{lo P 12 - fc) 2 } , 
= 2Dlp- 2 [u'{up' 2 - fc)]' + 4D 2 p~ 3 {iop" + lo'p'){lop 12 -k) + Dlp-\up' 2 - fc) 2 . 
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